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Fourier's Constants of Functions of Several Variables. 

By De. W. W. Kustekmann. 



In Vol. LVII of the MathematiscJie Annalen, A. Hurwitz has shown how 
to express the product of two ordinary Fourier's series in form of another 
Fourier's series, or, stated differently, how to compute Fourier's constants of 
the product of two functions from those of the functions themselves. In the 
following pages we solve the corresponding problem for multiple Fourier's 
series. In order to simplify matters the work is carried through for double 
Fourier's series only, but results and proofs admit of immediate and obvious 
generalization to «-tuple series. 

The solution of our problem depends vitally upon a proof of the relation* 



1 r^ /* T ™ 

71 •/— 7r«^— 7T /l,,P=0 



(a 2 + & 2 + c 2 + ^ 2 ), 



(1) 



where (a, b, c, d) A „ are Fourier's constants of f(x, y) and have the values 



=4 f f T f(«,p) 

71 *J —TC *J -TC 



II.V 



cos fj.a cosr/? 
sin j«a cos i'/? 
cos^a smv(3 
sin (ia sin »-/3 



dadfi. 



(2) 



The analogue of this relation for functions of a single variable is due to 
Parseval and was proved by him on the supposition that the function in ques- 
tion is really represented by its Fourier's development and that this series is 
integrable term by term. In 1893 de la Vallee-Poussinf gave a proof requiring 
merely that the function and its square be integrable. Hurwitz, in 1903, again 
called attention to the wide scope of the theorem, gave a new proof, and applied 

* In the following (a 2 +& 2 -j-e 2 +<J 2 )yai' stands for (a 2 +&V»+ C V''+^V'') and so on. E(z) denotes 

the largest positive integer contained in z. Hereafter we will write A/ip in place of 2 \]i+i) + \v+l) , 

(i when J u = i' = 0, 
2 " /x — O, x>0 or ^>0, x = 0, 
1 " ^>0, ^>0. 
•j- De la Vallee-Poussin, Ann. de la Soc.Scient. de Brnmelles, Vol. XVII, p. 18. 
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114 Kustermann : Fourier's Constants of Functions of Several Variables. 

it to the problem mentioned. More recently proofs of varying generality have 
been added by Fischer, Stekloff, Fatou,* and the formula has gained consider- 
ably in interest through the researches of Riess and Fischer (Riess-Fischer 
Theorem).! 

The first three paragraphs of our paper are devoted to proving relation (1) 
by a method which like those of de la Vallee-Poussin, Fischer and Fatou involves 
the use of Poisson's integral. (For a brief survey see §3.) We, too, make no 
assumption regarding the convergence of the Fourier's series of f(x, y). 

All we demand of f(x,y) is that it be bounded and possess a Rieinann 
double integral over the fundamental square.^ For convenience of language 
we will suppose f(x, y) periodic 2n, so that 



f(x±2nk, y±2nl)=f(x,y) for \ 



k=0,l,2,...., 
1=0,1,2, 



These conditions insure the existence of Fourier's constants (a,b,c, d)^, but 
do not in any way prejudice the question of convergence § of the formal 
Fourier's development for f(x,y) which we can write: 

00 

f{x,y)~ 2 l/A AJ ,[(a^cos^a;+&^ sm^.x)aosvy 

-\- (c^ cos (ix+d^ sin px) smvy]. (3) 

How relation (1) can be utilized to solve the problem originally proposed 
is shown in § 4. 

§ 1. Preliminary Theorems. 

Theorem I. If ${x,y) is bounded and integrable\\ in the fundamental 
square and can be represented by a trigonometric double series which is uni- 
formly convergent in this closed region, then the trigonometric series is a 
Fourier's series (that is its coefficients are given bv formulas of type (2) ). 

The proof is so closely analogous to the proof of the corresponding 
theorem for ordinary Fourier's series that we will not reproduce it here. 

* Fischer, Monatshefte, Vol. XV, p. 69; Stekloff, Mem. de I' Ac. de St. Petersbourg, Vol. XV, Ser. 8; 
Fatou, Acta Math., Vol. XXX. 

t Riess, "Uber orthogonale Funktionensysteme," Gott. Nachr., 1907; Fischer, Compt. Bend., 1907; 
also Young, Quart. Joum., Vol. XLIV; where this relation is shown to represent the necessary and 
sufficient condition that a given trigonometric series be the Fourier's series of a function whose square 
is integrable. 

t By the "fundamental square" we mean the square ( — it, — jr; v, ir) whose lower left-hand and 
upper right-hand vertices are ( — ir, — 7r) and (ir, ir), respectively. 

§ In the Stolz-Pringsheim sense. Of. Munchner Siteungsber., Vol. XXVII, p. 101. 

|| All integrals in this paper are to be taken in the Riemann sense. 
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Theorem II. If f(x,y) and ty{x,y) are both bounded and integrable in 
the fundamental square, if furthermore ty(x,y) can be developed into a double 
Fourier's series, uniformly convergent in this region, then 

~2 ) I f(v,y)<p(%,y)dxdy= 2 (aa+b(3 + cy + dS) liy / A^,* (4) 

where (a, b, c, d) ^ , (a, /3, y, 8) A „ are the Fourier's constants of f and <p, 
respectively. 

To prove this theorem multiply the Fourier's expansion of <p(x,y) by 
f{x,y) and integrate term by term over the fundamental square, f If, in 
particular, we let $(x, y) = f(x, y), equation (4) goes over into (1). The 
latter relation therefore holds in the special case, where / (x, y) can be 
developed into a uniformly convergent Fourier's series. To prove it for any 
integrable f(x,y) we need further preparations. 

Lemma: The double series 

S = 2 (a 2 + & 2 + c 2 + ^)„AV„, (5) 

where (a, b, c, d)^ are the Fourier's constants of a bounded integrable function 
f(x, y), is absolutely convergent and can be summed by rows, columns or 
diagonals. 

Peooe: Define ^(x,y) by the identity 

f ann 1 m 

^(x, y)=f(x, y)— | — + - 2 (a^ cos (ix+b^ sin^a) 

1 n 

+ - 2 (a 0v cosvy + c 0v sinvy) 

m n 

+ 22 [(a, Y cosfi'x+bp,,?, sin^'rc) oosv'y 



+ (c„,> p , cos (a' x + dp „, sin^'*) sin 1/3/] k 



* Cf. note * p. 113. 

t The reduction is effected by means of the well-known formulas: 



/"" sin sin, , f 0, if k:±\, 
/ Biniuo • cob \x die = 0; / *™ K! » da = 0, if k :£ 0. 
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Squaring both sides of this equation 



l^( x ,y)] 2 =[f( x ,y)T-\^f{x,y)+X(a li0 f(x,y)cosfix + b /i0 f(x,y)sinfix) 

A 11=1 



+ 2 (a „f(x,y) cos vy + c^ f(x,y) sinvy) 

v=l 
m n 

+ 2 2 2 [a^^ f(x,y) cos ^.'x cos v'y + b ll/ , v , f{x,y) sin ^.'x- cos v'y 



+ c^ v , f(x,y) cos fi'x sinv'y + d ll , v ,f(x, y) sin^'# sinv'y] 

1 m 

00 I X V (rfi 



/7 2 1 » 1 



(6) 



+ 77 + j 2 (a^cos^+fo^sin 2 ^) + - 2 (a 2 ^cos 2 ^/+c^sin 2 ^) 

m n 

+ 22 [a 2 -,- cos 2 ^'a3cos 2 T , '^ + & 2 -^ sink's; cos 2 v'y 

-\-c\, v , cos 2 [i'x sm 1 v' y -{- d\, v , sink's; sin* v'y] +2 J w2v2^'2»''. 
Here 2/« 2^ 2|ii' 2f' contains those product terms (except for constant coeffi- 
cients) which are of one of the sixteen types . ax • . vti • . u'x • c ? s v'y 
' ; ^ sm^ sin J sm ^ sin y 

where some of the numbers [i, v, u' , v' may be zero, but in no term of which the 
equalities u—u', v—v' coexist. If we now integrate equation (6) over the funda- 
mental square the terms resulting from 2/u 2v 2,u' 2v' vanish* and we have 

J f H(x,y)Tdxdy= j j [f(x, y) Ydxdy-rfl^ + 2 (a 2 +& 2 )„ 

n m n 

+ 2 (a 2 +c 2 ) ,+2 2 2 (a 2 +& 2 +c 2 +d S! )„„ 



a 2 I m 1 « 



+ t 2 |f + 2-S 1 (« 2 +^),o + r}2 i (a 2 + C 2 ) c 
+ 22 (a 2 +& 2 +c 2 +d 2 )4 

,14 = 1 J< = 1 J 

or 

-2 r r [^(^«/)] 2 ^% 

= ? 1 [/(^?/)] 2 ^^- 2 (a 2 + & 2 + c 2 + ^ 2 )„,/A /( ,. 

The integral on the left-hand side can clearly not be negative, hence, setting 

fi=m, v=n 



we can write 



8 mn = 2 (a 2 +& 2 + C 2 + i 2 )„,/A„ 
S mn <^-z f f [f(x,y)] 2 dxdy, 



* Cf. note f p. 115. 
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no matter what positive integral values in and n may have. Since, moreover, 
all the terms of (5) are positive, it follows from a well-known theorem* that 
this series converges absolutely and hence can be summed in the various ways 
mentioned. 

Cokollaby: The Fourier's constants of a bounded and integrable function 
of two variables approach zero when either one or both subscripts become 
infinite.^ 

§ 2. Discussion of Poisson's Double Integral. 

If in the double Fourier's series (3 ) we introduce the convergence factors 
r^ +v , where 0^r<l, we obtain a new double series 

00 iffr + V 

u r (x,y) = 2 — — [ (a ^ cos jix+b^ sin (/.x) cos vy 

+ (c^ cos [ix+d^ sin^a;) sin vy~\ (7) 

which is uniformly (and absolutely) convergent for all values of (x,y) and 
hence defines a function u r (x,y) which is finite, continuous, and integrable in 
the fundamental square. If we replace (a, b, c^d)^ by their values from (2), 
interchange summation and integration signs, and effect the summation by 
means of the formula 

|+j i ^cos^-|- 1 + r2 1 _7j cog<?) = ^(4»), 



we obtain 

u 



which by the transformation a — cc=%, (3 — 2/=>7 goes over into 
« r (x,y) = -L C Cf(x+£,y+ V )P r (£)PrWd£dy, 

= ttf rF(x,y,£,*i)P r (Z)PrWd£dn, (8) 

t 71 *S0 J 

where 

F{x,y,k,n)=f{xH, y+v)+f(x+£, y—n) +j '(x-£, y+ri)+f(x-£, y-n). (9) 

* Cf . Pringsheim, loo. cit., § 2, V and § 3. 

f This property of double Fourier's constants was proved from the integrals denning them by 
W. H. Young, Lond. Math. Soc. Proc, Ser. 2, Vol. XI, p. 133. 
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We wish to prove 

Theorem III. // f(x,y) is bounded and integrable in the fundamental 

square, then 

Urn u r (x,y)=f(x,y) 



r=l-0 



72) I 1 + I» + I S [. 



(10) 



at every point of continuity of f(x, y). 

Since 

c* /14-r s\ 

I P r {<p)d$ = 2arctan(— - — -tan-) <;?t 

and 

J-.7T 
P r (q>) d^> — 7t for every positive r<l, 


we may write 

u r (x,y)—f(x,y) 

= A f f \.F(x,y,£,y 1 )-4f(x,y)]P r tf)Pr(n)dZd>, 

1 f /*° /* a /* "" /*"' /^ a /'"' 1 

~4t1 2 \Jo Jo A ^0 Jo J« J 

~4ti 2 

We shall see that the sum of these three integrals can be made as small as we 
please by choosing r close enough to unity. 

\h\£f o a f Q a \F(v,y,Z,n)-±f(3>,y)\-\Pr(Z)\-\PrW\dZdri. 

Now 

1 — r 2 > 1— r 2 _ 1— r 

^rW)- 1 + r 2-2rcos<p= (1+r) 2 1 + r ^ 
for 0<[r<l; hence, |P r (<?>) | = P,($) ; furthermore, 

1 + r 2 — 2rcos$ = (l— r) 2 +4r sin 2 |>4r sin 2 -, 
for 0<r<l and 0<a<L$<L7t; therefore, 

P r ($>)<(l-r 2 )/4rsin 2 |. 

If (x,y) is a point of continuity for f(x,y), then a can be chosen so small 
that the upper bound of \F(x,y, £,yi) —<kf(x,y)\ in the domain Q = = a 

(11) 



is less than -\ . Thus, 



llil <«i for any <;»<!. 
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Denoting now by M the upper bound of | f{x,y) | in the fundamental square, 
we get 

|7,+/,|^8Jf| C P r (Z)d£ C P r {n)d n + CP r (£)d£ f P r {y,)dy,\ 

*> a ^0 ♦A) J a 

<4 7tM(l— r 2 )/rsin 2 ^, (12) 

Li 

a having been fixed so as to insure (11) we may choose jR<1 so close to unity 
that the right-hand side of (12) is less than e 2 for R^r<l; hence, 

\Ii\ + \h+h\<e 1 +e 2 , for R^r<l. 

It follows from (10) that lim u r (x,y)=f(x,y). 

r=l-0 

Corollaet: If the bounded function f(x,y) is integrable {according to 
Riemann's definition) in the fundamental square and u r (x,y) is defined by 
series (7), then 

lim -j f f [u r (x,y)] 2 dxdy=-j f C [f(x,y)] 2 dxdy. (13) 

Proof: Since f(x,y) is integrable in the sense of Riemann its points 
of discontinuity form a set of measure zero.* Consequently, [u r (x, y)] 2 
approaches [f(x, y)] 2 "almost everywhere."! Moreover, u r (x,y) is easily 
shown to be bounded in the aggregate for which <, r < 1. $ Hence, we can 
apply a theorem of Lebesgue § which states that if a sequence of integrable 
functions, bounded in absolute value, converges almost everywhere to a limiting 
function which is integrable, then the limit of the corresponding sequence of 
integrals is equal to the integral of the limiting function. It follows that (13) 
is true if the integrals involved are interpreted as Lebesgue integrals; but 
since, in view of our hypothesis, all these integrals have a meaning also 
according to Riemann's definition, and since the two kinds of integrals are 
identical, whenever they both exist, our corollary is proved for Riemann 
integrals. || 

* Cf. Hobson, "Theory of Functions of a Real Variable," § 311, p. 419. 
f /. e., except at a set of measure zero. 

t \Ur(w,y)\<^-4M ■ I P r (£)d£ I P r ( v )d v <M for 0<r<l. 

§ Lebesgue, Ann. de I'Ecole Normale, Ser. 3, Vol. XXVII (1910), p. 375. 

|| This method is no longer applicable when f(m,y) is only known to be integrable in Lebesgue's 
sense, for then the points of discontinuity may form a set of measure greater than zero. 



120 Kusteemann : Fourier's Constants of Functions of Several Variables. 

§ 3. Proof of the Generalised Parseval's Formula. 

We are now in a position to prove relation (1). Since the trigonometric 
double series 



r 



H + V 



u r (x,y)^ 2 —— [ (a ^ cos (ix+b^ sin fix) cos 



n, p=o A 



vy 



liV 



+ (c^ cosmic + d lip sin fix) sin vy J (7) 



is uniformly convergent we see from Theorem I, § 1, that it is the Fourier's 
development of u r (x,y) and that the Fourier's constants of u r (x,y) are 
r"+"(a, b,c,d)^. Hence, Theorem II is applicable and we have 

4 f f [u r (x,y)] 2 dxdy= 2 r ^ + "\a 2 +b 2 +c 2 +d 2 )J^ v . (14) 

7t */ — 7r V —IT fl, p=Q 

Now by the corollary of § 2, 

^f f U {x, y)T dxdy =\im -^ ( C [u r (x, y)] 2 dxdy 



= lim 2 r^+^ + ^ + c' + ^/A^. (15) 

r=l— #, j»=0 



The latter series of positive terms converges uniformly in r for the values 
0<r<gl, since it is term by term less than or equal to the series of positive 
constant terms (5) known to be convergent by the lemma of § 1. Hence, 
we may proceed to the limit r = l— in each term of (15) separately which 
gives (1). 

§ 4. Fourier's Constants of the Product f(x,y)-<p(x,y). 

It remains to be shown how relation (1) may be used to solve the problem 
originally proposed, to compute the Fourier's constants (21, $, ©, S))„„ of / • $, 
having given those of / and <p, which are (a, b, c, d)^ and (a, /3, y, 8)„„ 
respectively. Since / • <J> = i[/+<?>] 2 — 1[/— $Y and the Fourier's constants 
of /±$ are (a±a, b±(3, c±y, d±%) /iV we obtain from (1) 

^W^ I f-^>dxdy= 2 (aa+^+cy+dS^/A^. (16) 

To find 2I mre = — ^ I I / • ty cos mx cos w# (Za;^ we place for the moment 

71 J-icJ-t 

4 1 ' = $ cos ma? cos n?/. 
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The Fourier's constants of ty, namely (A', B', C, £>')„„, can be expressed 
in terms of those of $. For instance, 

A' llv = ~ g- ( ) ^'cos/iajcosr^/dajrfi/^^— j I ( 4>cosm:rcosp:rcosm/cosvi/da;<% 

1 Z* 71 " / ,1r 

— T - ! I I ^[cos(ff— ra)a; + cos( i M + w)a;] [cos(i>— «)«/+cos(v + w)2/]^ a:; ^ 

— 7\ G-ii—m, v—n ~T~ &ii,—m, v+n I ^+»», i>— » T~ G^t+m, v+n ) • 

Similarly, 

& liv — 4 \P/i— m, v—niP/i—rn, v+n T~ P/L+m, v—nrP/i+m, v+n) J 
l^/iv — 4 Kfii—m, v— n i /V— m, »+» i / /it+m, v—n~i /ii+m, v+n) ) 
"/«' — 4\')i-m, v—n I ® n—m, v+n \ ® n+m, v—n ~T ji+«i, v+b/ • 

But i|/' satisfies all conditions under which (16) was derived; hence, 
2U=Af f f'f-Vdxdy= I M' + &fl' + cC' + <M>'),>/A„. 

71 «^-ir J— ir n, v=0 

If we apply the same line of reasoning to 

4<" = $ sin mx cos «y, 



we find 



where 



Similarly, 
with 



and 



23 m „= v { a A"+bB"+cC"+dD") ILV /h ILV , 

in, v=0 

^liv — TV P/i—m, v—n P/L—m, v+nvPn + m, v—n~\P/i+m, v+n) > 
" [lv — i \^it— m, v—n~T~&ii.—m, v+n Q-n+m, v—n M/i+m, v+n) ) 
Is ILv — 4 \ Oft—m, v—n "#— m, »+n - r®/t+m, x— «"T" ®n+m, v+n) ) 
^/ik — 4 Vvjf— ro, »— n~r/ n— m, v+n //i+m, i»— ft //i+m, v+n) • 

e mtt = I (a.4'" + kB'"+cG v "+dZ)'")„/A„ 

A, l»=0 

-"/!>» — 4" V ^/it— Mj niTy/i-m, »+n / n+m, v—ni~/ii+m, v+n) > 
*J fi V — ~f \ ®tk—m, v—n \ Qft—m, v + n ®ii+m, v—n T" ^ii+m, v+n) j 
^/ij> ~~ 4 V^/t—wtj r— n dft—ni, v+n~T ^n+m, v—n ^fi+m, v+n) 5 
-^'/ttx — 4 \Pit.—m, v—n P/i—m, v + n i P/i+m, v—n P/L+m, v+n) } 



ft, »=0 



16 
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with 

■"■/lv — i V'/i- m, v— n ®/i— m, v+n O/L+m, v— n\~ Q/i+m, v+n) > 
■"pp — £ \ //»— m, v— »T/ h— m, v+n if n+m, v— n Yn + m, p,+n) > 
^ tiv 4\ P li— m, v— n I Pit,— m, v+n\ P n+m, v— n P/L+m, v+n) ) 

■L'/iv 4 \&n-m, v— n "ft— m, v+n &n + m. v— n I ^ii+m, v + n ) • 

Before using these formulas Fourier's constants with negative subscripts 
must be transformed into such with positive ones by considering that 

(a, /3, y, 3)_ K , X = (a, — [3, y, — %) kX , 

(a, &, y, S) _ K> _ x = (a, — /?, —y, B) kX , 

(a, (3, y, 5) „, _ x = (a, /?, — y, — 8) kX . 



